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Abstract 

We consider a single anharmonic oscillator with frequency u and coupling constant À respec- 
tively, in the strong-coupling regime. We are assuming that the system is in thermal equilibrium 
with a reservoir at temperature f3~ l . Using the strong-coupling perturbative expansion, we obtain 
the partition function for the oscillator in the regime À >> u, up to the order -t=. To obtain 
this result, we use of a combination of Klauder's independent- value generating functional (Acta 
Phys. Austr. 41, 237 (1975)), and the generalized zeta-function method. The free energy and the 
mean energy, up to the order -^=, are also presented. We are showing that the thermodynamics 
quantities are nonanalytic in the coupling constant. 
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1 Introduction 

The strong-coupling regime in quantum field theory is one of the unsolved problems of theoret- 
ical physics of the last century. There are many situations where one has to account for non- 
perturbative coupling regions, and to discuss the physics of the strongly coupled systems. Let us 
quickly mention a few examples. 

The standard case is QCD, where we believe there exists a deconfmement phase transition. It 
is expected that at sufficiently high temperatures, quarks and gluons become no longer bounded 
inside hadrons. There is also another phase transition which depends on the number of quark 
flavour: the chiral symmetry phase transition. The chiral symmetry is spontaneoulsly broken at 
zero temperature and it is restored at some finite temperature. Although QCD is asymptotically 
free in the ultraviolet limit, it has a coupling ~ 1 near the phase transition, where is the 

temperature. Consequently a non-perturative approach must be used to deal with this situation. 
The most fruitful idea is lattice QCD which is a means to study the phase transition from first 
principies. 

In condensed matter physics we also have different situations in which the strong-coupling 
regime may occur. For example, in the theory of superconductors there is the strong-coupling 
regime in the electron-phonon interaction. In particular, the strong-coupling limit for the polaron 
problem has been investigated by many authors. A detailed exposition can be found in Ref. 
[1]. Also in models with a disordered ground state one can use the strong-coupling approach. A 
standard case is the O(N) non-linear sigma model, describing lattice magnets, where there is a 
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regime in which one can treat the coupling between different sites as a perturbation. Until now we 
pointed out that there are many situations in quantum field theory that cannot be described using 
the weak-coupling perturbation theory. Nevertheless, the use of the perturbative expansion outside 
the weak-coupling regime is also fundamental for our understanding of the whole perturbative 
renormalization program in quantum field theory. 

Various types of harmònic oscillator have served for long as simple analogue systems for more 
general and complex situations in quantum field theory. Actually, the anharmonic oscillator, with 
a ji x 4 (t) term is formally very similar to the field theory describing a scalar field with a quartic 
self-interaction. Therefore, in this paper we present a method for calculating the partition function 
and the Helmholtz free energy for a single oscillator with the anharmonic ^ x 4 (r) contribution in 
the strong-coupling regime up to the order In other words, our model belongs to the regime 
in which the coupling constant is much larger than the frequency, i.e., À >> uj. This is also a con- 
venient way to study the regime of large fluctuations. In fact, the oscillator with any anharmonic 
term of the kind ^yy x 2p (t), p > 3 in the strong-coupling regime can be analysed by our method 
and also the generalization for N non-interacting anharmonic oscillators is straightforward. 

The anharmonic oscillator is described by the Hamiltonian 

and for simplicity we are assuming that our system is one-dimensional and is in thermal equilibrium 
with a reservo ir at temperature We are working in the imaginary time formalism and making 
use of the Kubo-Martin-Schwinger (KMS) condition [2] [3] [4]. 



To find the partition function and the Helmholtz free energy, our approach consists in the 
combination of two techniques used currently in the literature: the strong-coupling expansion [5] 
[6] [7] [8] [9] [10] and the zeta-function method [11] [12], which is a classical case of a general 
formalism concerning regularized products and determinants. For a review of the strong-coupling 
expansion in field theory, see for example Ref. [13] and also Ref. [14], and for a recent treatment 
of the strong-coupling expansion in quantum mechanics, see for example Ref. [15]. 

Let us briefly discuss the strong-coupling expansion in Euclidean field theory at zero tempera- 
ture. The bàsic idea of the approach is the following: in a formal representat ion for the generating 
functional of complete Schwinger functions of the theory Z(h), we treat the Gaussian part of the 
action as a perturbation with respect to the remaining terms of the functional integral, i.e., in the 
case for the (\ip A )d theory, the local self-interacting part, in the functional integral. In the litera- 
ture this approach it has been called the strong-coupling expansion. A stimulating reading can be 
found in Ref. [16] and Ref. [17]. The main difference from the the Standard perturbative expan- 
sion is that we have an expansion of the generating functional of complete Schwinger functions in 
inverse powers of the coupling constant. We are developing our perturbative expansion around the 
independent- value generating functional Qo(h), where different points of the Euclidean space are 
decoupled, since the gradient terms are dropped [18] [19] [20] [21]. For study the analytic structure 
of the independent-value generating-functional in the complex plane of the coupling constant À, 
see for example Refs. [22] and [23] . 

In the perturbative expansion around the independent-value model, the representation of the 
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Schwinger functional Z(h) is defined by an unrenormalized perturbative series, which can be 
truncated in the order of the approximation. For example, if Z(h) = Qo(h), we call it the 
independent-value approximation or the zeroth-order approximation. It should be noted that 
although the independent-value functional Qo(h) is not a product of Gaussian integrals, it can 
be viewed formally as an infinite product of ordinary integrals, one for each point of the d- 
dimensional Euclidean space. The fundamental problem of the strong-coupling expansion is how to 
give meaning to the independent-value generating functional and to the unorthodox representation 
for the Schwinger functional. 

A solution to this problem, was presented by Klauder and colaborators a long time ago [19] [24] 
[25]. They obtained a quite interesting expression for the independent-value generating functional 
describing Euclidean free and also self-interacting scalar fields. We would like to stress that a 
naive use of a continuum limit of the lattice regularization for the independent-value generating 
functional leads to a Gaussian theory, where we simple make use of the central limit theorem. 
The fundamental modification which allow us to avoid the central limit theorem is a change in 
the measure in the functional integral. Consequently, in the Klauder's derivation, starting from 
the independent-value generating functional, and taking the limit where the coupling constant 
becomes null, the interacting theory solutions become the pseudo-free one. 

In this paper we show how it is possible to compute the partition function and the Helmholtz 
free energy of the anharmonic oscillator in the strong-coupling regime, up to the order Our 
results are showing that the thermodynamics quantities are nonanalytic in the coupling constant. 
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We would like to stress that although, by a unitary transformation the quartic anharmonic oscil- 
lator can be expressed in a form where the dependence on the coupling constant appears in an 
overall factor, we are using the Standard form, since we are interested in the regime where the 
anharmonic terms are much larger than the kinetic term. In other words we are interested in the 
regime of large fluctuations. 

This paper is organized as follows: In section II, the strong-coupling expansion for a single 
anharmonic oscillator is presented. In section III we calculate the partition function and other 
thermodynamics quantities that are derived from it, as the Helmholtz free energy and the mean 
energy of the system. Finally, section IV contains our conclusions. We assume that the physical 
quantities are dimensionless. Consequently, it is convenient to introduce an arbitrary parameter 
H with mass dimension to define all dimensionless physical quantities. For simplicity, we assume 
that fi = 1 since we are not interested in the scaling behavior of the model. As usual (3' 1 is the 
temperature of the thermal bath. 

2 The partition function of a one-dimensional quantum 
mechanical system 

Let us consider a one-dimensional quantum mechanical system. The partition function for the 
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system assuming that it is in thermal equilibrium with a reservoir at temperature 1 is given by 



where in the functional integral we require that x(r) is periòdic with period (3, i.e., x(r) = x(t+(3). 
There are many different physical situations that can be analysed starting from the partition 
function. We would like to discuss the case of a single anharmonic oscillator, where the contribution 
of V(x(t)) is given by 



For simplicity, we are choosing m 2 — 1. As we will see, it is not difficult to apply our method in 
non-polynomials theories. The second situation that can also be analysed is the non-polynomial 
model, defined by the following expression 



In order to study the anharmonic oscillator, let us sketch the solution for the single harmònic 
oscillator. In the potential Vi(x), choosing À = we obtain the harmònic oscillator, where the 
partition function can be found in texbooks. The partition function can be written in a more 
tractable way, where an integration by parts is done and a surface term is disregarded. To 
avoid minor complications, we will assume the following expression to the free partition function 
Z(/3)\\=o. In particular, we shall therefore start from 



Z((3) = í [dx(r)\ exp f- f dr ( A 2 + V(x(r))) 

Jx(p)=x(p) Jo \2 dr j 



(2) 




(3) 




(4) 



Z([3)\ x=0 = í \dx(r)} exp - / dr -x(r) 

Jx(0)=x(3) JO 2 



( 



dr 2 



+ uj 2 x(r) 



(5) 
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It is a Standard procedure to define the following kernel K(u; r — t') by the equation 

/ d 2 



K(cu; r - t') 



dr 2 



(6) 



and substituting Eq.(6) in Eq.(5) the free partition function Z(f3) becomes 



Z({3)\ x=0 = f [dx(r)\ exp 

Jx(0)=x(/3) 



f/3 r/3 1 

/ dr dr' -x(t)K(u;t - t')x(t') 
Jo Jo 2 



(7) 



As usual, we define the generating functional Z{(5] h) introducing an external source h{r). At this 
point it is convenient to consider h(r) to be complex. Consequently h{j) = Re(h) + ilm(h). In 
the paper we are concerned with the case Re(h) = 0. Therefore the generating functional Z((3; h) 
is defined by 



Z(P;h)\ x=0 = í [cfe(r)]exp - F dr F ' dT J L·(T)K(u; t-t')x(t')+ FdTh(r)x(r) 

Jx(0)=x(fí) Jo Jo 2 Jo 



(8) 



Note that we are using the same notation for functionals and functions, for example Z(f3;h) 
instead of the usual notation Zp[h]. Since the integrations that appear in Eq.(8) are Gaussian it 
is straightforward to write 



Z(p; h)\x=o = Zp exp 



/ dr í dr' -h{r)G{uJ- r - r')h{r') 
Jo Jo 2 



(9) 



where the partition function is defined by Z(f3) = Z(f3; h])h=o, and the Green function G(u; r — r') 
is the inverse kernel, defined by 



F dr' K(u,t - r')G(u,r' - r") = 5(r - r"). 
Jo 
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At this point it is important to define the modified kernel K(u, a;r — r') by the equation 

K(u, a- r-r')= (-^ + (1 - a) co^j 5(r - r'), (11) 

where a is a complex parameter defined in the region < Re (a) < 1. Note that in general, 
Re (cr) 7^ 1, because Re(cr) = 1 introduce infrared divergences in the calculations. The zero 
frequency case (cr = 1) in the modified kernel can be assumed in some very special situations 
only, as for example to calculate the renormalized vacuum energy of a scalar field in the presence 
of boundaries, where Dirichlet boundary conditions are assumed. See for example Ref. [14]. 
To summarize, the choice of a suitable o will simplify our calculations in some situations. This 
modification will be clarified in the next section. 

Let us now suppose an anharmonic oscillator in the strong-coupling regime. We would like 
to stress that the semiclassical approximation (the WKB approximation) can also be used in 
the strong-coupling regime. The difliculty is to find the classical orbits which are stationary- 
phase points in the functional integral [26]. First, as we made in the non-interacting case, it is 
convenient to couple linearly the anharmonic oscillator to a r-dependent external source. Therefore 
the generating functional at finite temperature is given by 



Z(/3;h)= f Jdx(r)} 

Jx(ti)=x(/3) 

- f dr f 13 dr' -x(t) K x(t') + f dr { --au 2 x 2 {r) - -x\r) + h{r)x{ 
Jo Jo 2 Jo \ 2 4! 



where we have also integrated over all periòdic paths and K = K(üú,o~;t — r'). Functional 
differentiation gives the thermal average of a time-ordered of position operators, i.e., the correlation 
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functions for a stochastic process. For sake of completeness, we would like to present the simple 
result for the partition function of the anharmonic oscillator in the regime À << u>, in first order 
in À. One finds [27] 

Z(P) = (2smh(^)) _1 (l - | [ dri(coth 2 (^)) + 0(À 2 )) . (13) 

To find the partition function for the anharmonic oscillator in the strong-coupling regime it is 
natural to use an unorthodox perturbative theory, i.e., the strong-coupling perturbative expansion. 
The idea is to treat the Gaussian part of the action as a perturbation with respect to the non- 
Gaussian terms in the functional integral. We get the following formal representation for the 
generating functional at finite temperature Z(/3; h): 

z ®» = ^{-k£*Jo'"sm K{a · r > T -*sfó) w ·' T! ' i, ' (14) 

where Q((3,a; h), the new independent- value functional integral, is given by 

Q(/3, <x ; h) = JV í [dx(r)} exp f dr ( -\a uo 2 x 2 (r) - ^-x\t) + h(r)x(r) ) , (15) 

Jx(o)=x(p) [Jo \ 2 4! J 

and the modified kernel K(u, a;r — t') was defined by Eq.(ll). The factor J\í is a normalization 
that can be found using that Q{(5, a; h)\h=o = 1. 

The main difference from the Standard representation for the generating functional is that we 
have an expansion of the generating functional in inverse powers of the coupling constant. We 
are developing our perturbative expansion around the independent-value generating functional 
Q(f3,<j] h). We would like to stress that we are considering a modification of the strong-coupling 
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expansion. We split the quadratic part in the functional integral, which is proportional to the 
frequency squared, into two parts; one contributes together with the derivative term in the action 
as the perturbation, and the other appears in the independent-value generating functional. 

One way to proceed is to neglect high-orders terms in the perturbative expansion. Therefore 
in the leading order, we have that Z(f3;h) can be written as 

Z(ft h)=(l-l f ér [ dr-^K^a; r - T <)^) W, <r, h), (16) 

Since we are mainly interested in presenting the partition function, we can also assume that the 
external source is constant i.e. h{r) = h. As we will see, this assumption will lead us to redefine 
Klauder's representation for the independent-value generating functional. Of course, the above 
assumption considerably simplifies our problem, but we still have some work to do. 

To evaluate ln Z((3; h), note that we have two steps to follow. The first one is to give meaning to 
the independent-value generating-functional, and the second one is to regularize and renormalize 
the kernel K(u,a; r — t') integrated over the volume [0, /?]. Note that the parameter a was 
introduced only to simplify our calculations in some situations. Therefore a can be complex if we 
are able to work in all order of perturbation theory. The generating functional does not depends 
on the value for a. Since we concentrate in the leading order, some care has to be taken to prevent 
a complex generating functional. A simple way to avoid the problem is assume that the parameter 
o is real. Therefore we will impose that Im (cr) = 0. In this situation the independent-value 
functional Q(f3, cr; h) should be a normalized, positive defmite functional. We will discuss this 
point latter. 
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In the next section, we will use a combination of Klauder's independent-value generating 
functional, and also use the generalized zeta-function method to regularize and renormalize the 
kernel K(u, a;r — r') integrated over the Euclidean time. 

3 The partition function for the anharmonic oscillator in 
the strong-coupling regime 

In the present section we study the single anharmonic oscillator in the strong-coupling regime. 
Since all the thermodynamics quantities are derived from the Helmholtz free energy, let us proceed 
in deriving the free energy. The Helmholtz free energy can be obtained from ln Z(/3, h)\h=o, i.e. 
F{(3) = —jj\nZ(f3; h)\h=o- To have a well defined meaning to the Helmholtz free energy that can 
be obtained from Eq.(16) we may proceed as follows. 

First, we use the Klauder's result, as the formal definition of the independent-value generating 
functional Q(f3, a; h] . Second, we have to regularize and renormalize the kernel K(lü,<j;t — r') 
integrated over the Euclidean time. The sucess of our method, depends critically on the possibility 
to handle the independent-value generating functional and the kernel integrated on the Euclidean 
time [0,(3]. 

Since we are concerned with the strong-coupling regime, to evaluate InZ(fl) let us use the 
leading term. Using the cumulant expansion idea, which relates the mean of a exponential to the 
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exponential of means, after some simple calculations we obtain 



where £(s) is the global generalized zeta-function associated with the operat' 




dr 2 



+ (l-a) u 2 ). 



(17) 



There are some issues that we would like to discuss. Note also that, we assume triennal equi- 
librium and since we are working in the Euclidean formalism, the spectrum of the operator 
D = (— -£2 + (1 — <y)oJ 2 ^J has a denumerable contribution. Remind that o is a complex parameter 
defined in the region < Re (a) < 1. At this point, let us impose that Im (a) = 0, since Q(/3, a; h) 
should be a positive defmite functional, i.e., a characterized functional of a generalized stochastic 
process. 

The operator D is a positive defmite eliptic operator, and has a complete set of orthonormal 
eigenfunctions x n (r) and associated eigenvalues a n . We have 



with the boundary conditions x n (0) = x n (j3). Note that we have fí f dr x n (r)x n i(r) = 5 nn r. 

At this point, we would like to stress that a number of difïiculties appear for complex o. 
Since we are working in the leading order of the inverse power of the coupling constant, another 
problem related to the choice Im (<r) 7^ 0, is related to the eigenvalue equation given by Eq.(18). 
For complex a, the eigenvalue equation involve complex eigenvalues a n and the eigenfunctions still 
form a complete, but not orthogonal set. See for example [28] [29]. Actually there is a program to 
study non-self adjoint operators and this is related to the Parisi-Wu stochastic quantization [30] 




(18) 
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with a complex Langevin equation [31] [32]. Actually complex Euclidean action also appears in 
different systems. For example the Euclidean action for a gauge theory with external static charge 
is complex. Effective actions with topological terms are complex in the Euclidean formalism. 
Such kind of problems also appear in lattice QCD, where the analysis of the phase diagram in 
the temperature-chemical potential plane has been investigated [32]. It is well known that in the 
SU(3) theory, the fermionic determinant is complex, giving to a complex Euclidean action. For 
the study of the Gross- Neveu model with a nonzero imaginary chemical potential, see Ref. [33]. 
Summarizing, complex Euclidean generating functional introduce new stimulated problems and 
deserves a carefull analysis of them. 

Going back to the eigenvalue equation, the generalized zeta-function associated with the op- 
erator (--£; + (1 - a)uj 2 \ i.e., (_ ^ is defined by 



dr 



(19) 



where the spectrum is given by 



/27mV 



+ (1 -u)üj 7 



, ne Z. 



(20) 



Using the definition for the global generalized zeta-function and the spectrum of the operator 
given by Eq.(20) we have that the generalized zeta-function is given by 



C-JL +(i- CT)w 2( s ) - E 



2nn 

7T, 



;i - a) uj 2 



(21) 



Here, it is useful to define the modified Epstein zeta-function in the complex plane s, i.e., the 
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function Ç(s, v) by: 

oo 

C(8,u)= E (n 2 + v 2 r, v 2 >0. (22) 

n=— oo 

The series defined by Eq.(22) converges absolutely and defines in the complex s plane an analytic 
function for Re(s) > |. It is possible to analytically extend the modified Epstein zeta-function 
where the integral representation is vàlid for Re(s) < 1, [34] [35]: 

oo 
n=—oo 

For a different representation for the analytic extention of the modified Epstein zeta-function in 
terms of the modified Bessel function K a {z) or the Macdonald's function, the reader can use for 
example Ref. [36]. 

It is not difficult to write the generalized zeta-function in terms of the modified Epstein zeta- 
function. We have 

<-£ +(1 -,>>>=(|)^(s.^(^)). <*) 
where Ç(s, v), is the modified Epstein zeta-function. As we discussed, the series representation for 
Ç(s, v) converges for Re(z) > ^ and its analytic continuation defines a meromorphic function of s 
which is analytic at s = 0. The modified Epstein zeta-function has poles at s = |, — |, etc. Using 
Eq.(23) is not difficult to show that the vàlues for the modified Epstein zeta-function ((s,i/), at 
s = and -§^Ç(s, ^)| s =o are given by 

Ç(s,v)\ s=0 = 0, (25) 



+ 4 sin 7ts 



i: 



(t 2 



s dt 



(23) 
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and also 



d_ 

ds 



í/)| s=0 = — 21n(2 sinh7rz/). 



(26) 



Since we are interested in calculating the derivative of the generalized zeta-function at the origin 
of the complex s plane, we have 

2s 



ld_ 

2fc 



1 / p \ 2s d 



2 \ 2vr / ds 



CM 



(27) 



Choosing (T = 0, and using Eq.(25) and Eq.(26) in Eq.(27) we obtain the well known result in the 
literature. For the general case (cr ^ 0) we have 



--c d2 

2 ds 



+ (l-cr)w : 



s=0 



-ln 



(2 sinh ((l-*)^f 



(28) 



Thus it is clear that the zeta function regularization can be used to control the divergences of the 
kernel K(u>, cr; r — r') integrated over the Euclidean time. 

We would like to stress that we are using the Klauder's result, as the formal definition of 
the independent- value generating functional derived for scalar fields in a <i-dimensional Euclidean 
space. We would like to point out that in Klauder's derivation for the independent-value model a 
result was obtained which is well defined for all functions which are square integrable in R n i.e., 
h(x) e L 2 (R n ). This observation allow us to conclude that we need also to use a normalization in 
the situation that we are investigat ing. It is possible to show that the independent-value generating 
function can be written as 



Q((3, cr; h) = exp 



2(3 Jo 



fi 



dr 



00 du 



oo \U\ 



1 — cos(/wí)) exp 



1 2 2 ^4 
a üü U -U 

.2 4! , 



(29) 
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There is no need to go into details of this derivation. The reader can find it in Ref. [16] [19] [25]. 
It is important to stress that the independent- value generating function defined by Eq.(29) do not 
reduce to the conventional free independent- value model, if we choose À = 0. The non-gaussian 
contribution represents a discontinuous perturbation of the free theory, or using the Klauder's 
definition; a pseudo-free theory. Actually, this is the key point of the program developed by 
Klauder to investigate non-renormalizable models in field theory. For a interesting study of the 
pseudo-free harmònic oscillator see for example Ref. [37]. 

In order to study Q(j3, a; h), let us define E(u, a, A; h) given by 

/oo du í \ A \ 

1—7 (1 - cos(hu)) exp — cruj 2 u 2 - — u A ) . (30) 
-oo \u\ \ 2 4! ) 

Using a series representation for cos x, it is not difficult to show that 

E(u>, a, A; h) = 2 £) jT duu 2k ~ l exp (-1 ou? u* - A . (31) 

Now let use the fact that the a parameter can be choosen in such a way that the calculations 
becomes tractable. Analysing only the independent-value generating functional it is not possible 
to write Q(/3, a\ h) in a closed form even in the case of constant external source. One way to obtain 
a closed expression is to choose a = 0. Therefore we have 

00 ( — 1^ /-oo \ 

^^A^jU^^^I duu^eM—u 4 ). (32) 
At this point let us use the following integral representation for the Gamma function [38] 



f-oo | v I y\ 

\ dxx u - 1 exp(- l· ix p ) = -yTv Y - , Re(//) > Re(z^) > p > 0. (33) 
Jo V \VI 
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It is clear that the (Àx p ) theory, for even p > 4, can also easily handle applying our method. 
Using the result given by Eq.(33) in Eq.(32) we have 

E(u;,a,\;h)\ a=Q = Y,9(k)-T, (34) 

k=l A2 

where the coeficients g(k) are given by 

**> = 5Sr (4l),r( 5 ) · (35) 

Substituting the Eq.(34) and Eq.(35) in Eq.(29) we obtain that the independent- value generating 
function Q(P, a\ h)\ a=0 can be written as 

r p 00 tfk' 



Q(P,<t; h)\ a=0 = exp 



1 rP 00 h AK 

2p Jo ^ A2 



(36) 



It is easy to calculate the second derivative for the independent- value generating function with 
respect to h. Note that Q(/3,a; h)\ h=a=0 = 1. Thus we have 

^ U=o = 2^ 9{k){2k)(2k - J ex P L· J + ( 37 ) 

where G(h) is given by 

(00 L2fc+2<j-2\ / 1 00 h 2k \ 

E Í/(M) — ex P lE#)Tr . ( 38 ) 
fc,,=l A 2 J \ l k=1 \2 ) 

and g(k,q) = kqg(k)g(q). We are interested in the case h — 0, therefore the double series does 
not contributes to the Eq.(37), since limh^oG(h) = 0. Using the fact that we are interested in the 
case h = 0, we have the simple result that in the Eq.(37) only the term k — 1 contributes. We get 

d 2 Q((3, a\ h) /3tt 



dh 2 lh - a -° V 8A' 
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(39) 



Substituting the result obtained from the generalized zeta-function method given by Eq.(28) 
(choosing cr = 0) and Eq.(39) in Eq.(17) we have that lnZ(/3) is given by 

(40) 



i-ln( 2s mh(íf) 



Therefore the partition function for the single oscillator is 



1 /3F 
e2 V 8A 



(2sinh(f); 

It should be noted that for a system of N harmònic oscillators, the partition function is 



(41) 



N ■ 



(42) 



(2sinh(f ; 

Therefore the identification N = y^ff , the partition function for the strongly coupled single 
oscillator seem to be proportional to the partition function of N harmònic oscillators, up to the 
order 

Other thermodynamics quantities that we are able to find are the Helmholtz free energy and 
the mean energy. The Helmholtz free energy is given by F(f3) = — jj \nZ(/3, h)\h=o- Thus we have 



1 3tt 



-- + ln 2smh(^-; 



P V 8A 

It is possible to write this expression in the following way 



(43) 



(44) 



Finally the mean energy is defined by E = — ^ ln Z((3)\ h=0 . Therefore we have 



E 




LO 



+ 



LO 



(45) 
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We conclude with some observations. Usually, the free energy of a finite system is an analytic 
function of the parameters that define our physical system. Nevertheless, our results are showing 
that the thermodynamics quantities are nonanalytic in the coupling constant. 

The picture emerging from the previous discussion is the following: in the strong-coupling 
perturbative expansion we may split the problem of defining the generating functional into two 
parts: how to define precisely the independent-value generating functional and how to go beyond 
the independent-value approximation, taking into account the perturbation part. Our results show 
that the strong-coupling perturbative expansion, in combination with an analytic regularization 
procedure, is a useful method to compute global quantities, as the Helmholtz free energy, in the 
strong-coupling regime. 

4 Conclusions 

In this article we studied the strong-coupling regime in one-dimensional models, after analytic 
continuation to imaginary time. One-dimensional models are very simple system for which we 
can apply our method in obtaining thermodynamics quantities in the leading order in the inverse 
of coupling constant. We calculate the partition function and the Helmholtz free energy for the 
anharmonic oscillator, using the strong-coupling perturbative expansion and the generalized zeta- 
function analytic regularization. It was possible to present expressions up to the order for the 
partition function and the other thermodynamic quantities derived from the the Helmholtz free 
energy. 
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The picture that emerges from our method is the following: in the strong-coupling perturbative 
expansion we may split the problem of defining the generating functional into two parts: the flrst 
is how to define precisely the independent-value generating functional. We may use of a lattice 
approximation to give a mathematical meaning to the non-Gaussian functional. Actually, it is not 
easy to recover the continuum limit, and as we discussed the use of the central limit theorem leads 
us to a Gaussian theory. Instead of this, we are using the Klauder's independent-value generating 
functional. This is the key point of the article and can explain the discrepance between our results 
and the Ref. [27] (pp. 930, 931), where the large coupling behavior of the anharmonic oscillator is 
discussed. 

The second part is to go beyond the independent-value approximation and take into account 
the perturbation part. This problem can be controled using an analytic regularization. Besides 
these technical problems, we still have the problem of obtaining the Green's functions from this 
approach. The strong-coupling perturbative expansion is not fit to obtain local quantities, as the 
Green's functions of the model. On the other hand, our results show that the strong-coupling 
perturbative expansion, in combination with an analytic regularization procedure, is a useful 
method to compute global quantities, as the Helmholtz free energy, in the strong-coupling regime. 

There are several directions for investigations, using the Klauder's result and an analytic 
regularization procedure. It should be possible for applying the method to more realistic theories. 
To mention a few: since scalar fields play a fundamental role in the Standard model, the study of 
the strongly coupled (A</? 4 )d theory at finite temperature and also the renormalized vacuum energy 
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of scalar quantum fields in the presence of macroscopic structures deserves future investigations. 
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